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Bases of finite-dimensional Hilbert spaces (in dimension d) of relevance for quan- 
tum information and quantum computation are constructed from angular momentum 
theory and su(2) Lie algebraic methods. We report on a formula for deriving in one 
step the (1 + p)p qupits (i.e., qudits with d = p a prime integer) of a complete set 
of 1 + p mutually unbiased bases in C p . Repeated application of the formula can 
be used for generating mutually unbiased bases in with d = p e (e > 2) a power 
of a prime integer. A connection between mutually unbiased bases and the unitary 
group SU(d) is briefly discussed in the case d = p e . 
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I. PROLEGOMENA 

The use of symmetry adapted functions (or state vectors) is of paramount importance in 
molecular physics and condensed matter physics as well as in the clustering phenomenon of 
nuclei. For instance, wavefunctions adapted to a finite subgroup of SU(2) turn out to be 

n n 

very useful in crystal- and ligand-field theory 

It is the purpose of the present work to report on state vectors adapted to the cyclic 
subgroup Cd of 5*0(3) ~ ST/(2)/Z 2 . Such vectors give rise to bases of SU(2), the so-called 
mutually unbiased bases (MUBs), to be defined in Section V. They play a fundamental role in 
quantum information and quantum computation in view of the fact that these bases describe 
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qudits (the analogs of qubits in dimension d). A single formula for MUBs is obtained in this 
paper from a polar decomposition of SU{2) and analysed in terms of quantum quadratic 
discrete Fourier transforms, Hadamard matrices, generalized quadratic Gauss sums, and 
special unitary groups SU(d). 

Before dealing with the main body of this work, we continue this introduction with a brief 
survey for nonspecialists of some particular aspects of quantum computation and quantum 
information for which the concept of MUBs is useful. 

According to the law by Moore, the size of electronic and spintronic devices for a classical 
computer should approach 10 nm in 2018-2020, i.e., the scale where quantum effects are 
visible, a fact in favor of a quantum computer. This explains the growing interest for a new 
field, viz., the field of quantum information and quantum computation. Such a field, which 
started in the 1980 's, is at the crossroads of quantum mechanics, discrete mathematics and 
informatics with the aim of building a quantum computer. We note in passing that, even 
in the case where the aim would not be reached, physics, mathematics, informatics and 
engineering will greatly benefit from the enormous amount of works along this line. 

In a quantum computer, classical bits (0 and 1) are replaced by quantum bits or qubits 
(that interpolate in some sense between and 1). A qubit is a vector \<f>) in the two- 
dimensional Hilbert space C 2 : 

\(/>)=x\0)+y\l), xeC, yeC, \x\ 2 + \y\ 2 = 1, (1) 

where |0) and |1) are the elements of an orthonormal basis in this space. The result of a 
measurement of |0) is not deterministic since it gives |0) or |1) with the probability \x\ 2 or 
\y\ 2 , respectively. The consideration of N qubits leads to work in the 2 Af -dimensional Hilbert 
space C 2JV . Note that the notion of qubit, corresponding to C 2 , is a particular case of the 
one of qudit, corresponding to C d (d not necessarily in the form 2^). A system of N qudits 
is associated with the Hilbert space C . In this connection, the techniques developed for 
finite-dimensional Hilbert spaces are of paramount importance in quantum computation and 
quantum calculation. 

From a formal point of view, a quantum computer can be considered as a set of qubits, 
the state of which can be (controlled and) manipulated via unitary transformations. These 
transformations correspond to the product of elementary unitary operators called quantum 
gates acting on one or two qubits. Measurement of the qubits outcoming from a circuit of 
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quantum gates yields the result of a (quantum) computation. In other words, a realization 
of quantum information processing can be performed by preparing a quantum system in a 
quantum state, then submitting this state to a unitary transformation and, finally, reading 
the outcome from a measurement. 

Unitary operator bases in C d are of pivotal importance for quantum information and 
quantum computation as well as for quantum mechanics in general. The interest for unitary 



operator bases started with the seminal work by Schwinger [121 ] . Among such bases, MUBs 
play a key role in quantum information and quantum computation 13j- 17J. From a very 
general point of view, MUBs are closely connected to the principle of complementarity in- 
troduced by Bohr in the early days of quantum mechanics. This principle, quite familiar in 
terms of observables like position and momentum, tells that for two noncommuting observ- 
ables, if we have a complete knowledge of one observable, then we have a total uncertainty 
of the other. Equation ( BTj) in Section V for a ^ b indicates that the development in the 
basis B a of any vector of the basis Bj, is such that each vector of B a appears in the devel- 
opment with the probability 1/d. This is especially interesting when translated in terms of 
measurements, the bases B a and B\> corresponding to the (nondegenerate) eigenvectors of 
two noncommuting observables. 



II. A NONSTANDARD BASIS FOR SU(2) 

The various irreducible representation classes of the group SU (2) are characterized by a 
label j with 2j G N. The standard irreducible matrix representation associated with j is 
spanned by the irreducible tensorial set 

B 2 j+i = {\j,m} ■ m = j, j - 1, . . . , -j}, (2) 

where the vector \ j, m) is a common eigenvector of the Casimir operator j 2 and of the Cartan 
operator j z of the Lie algebra su{2) of SU(2). More precisely, we have the relations 

j 2 \j,m) = j(j + 1)|.7, m), j z \j,m) = m\j,m), (3) 

which are familiar in angular momentum theory. (We use lower case letters for operators 
and capital letters for matrices so that j 2 in ([3]) stands for the square of a generalized angular 
momentum. ) 
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Following the works in 18|, [19( , let us define the linear operators v ra and h by 



e i2 ^|j,-j)(j,j|+ E g°'- w)a |i,m + l)(i, 



m 



(4) 



m=-j 



and 



/i= E vO' + m )0'- ? " + 1 )U>" i )0'> m l> 



(5) 



m=-j 



where 



r G R, q = e 



2l?i/(2j+l) 



a G R. 



It can be checked that the three operators 



1 r 



j+ = hv ra , j'_ = (v ra ) j h, j z = - \h 2 - (v ra yh 2 v, 
where (v ra y stands for the adjoint of v ra , satisfy the commutation relations 

[jz , j+] = +j+ , [jz ,j-] = -j- , [j+ ,j-] = 2jz 



(6) 



(7) 



(8) 



of the algebra su(2). (In angular momentum theory, the operators j + and j_ are connected 
to j 2 via j 2 = j±j T + j z (j z =F 1).) 

The operator v ra is unitary while the operator h is Hermitian. Thus, Eq. (jTJ) corresponds 
to a polar decomposition of su(2) with the help of the operators v ra and h. It should be 
noted that v ra can be derived in terms of operators acting on the tensor product of two 
commuting quon algebras associated with two truncated harmonic oscillators. The latter 



k- 
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ermions which are supersymmetric 



oscillators play a central role in the intoduction o: 
objects interpolating between fermions and bosons 

It is evident that v ra and j 2 commute. Therefore, the {j 2 ,v ra } scheme constitutes an 
alternative to the {j 2 ,j z } scheme. This yields the following result. 



Result 1. For fixed j , r and a, the 2j + 1 vectors 

1 j 



\ja; ra) 



1^ 

m=-j 



(j+m)(j—m+l)a/2—jmr-\-(j+m)a 



(9) 



with a = 0,1,..., 2j, are common eigenvectors of v ra and j 2 . The eigenvalues of v ra are 
given by 



v ra \ja;ra) = q^ a+r ^ a \ja;ra), 



(10) 
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so that the spectrum of v ra is nondegenerate. 

For fixed j, r and a, the inner product 

(ja;ra\j(3;ra) = 5 a: (3 (11) 

shows that {\ja; ra) : a = 0, 1, . . . , 2j} is an orthonormal set which provides a nonstandard 
basis for the irreducible matrix representation of 577(2) associated with j. 



III. QUANTUM QUADRATIC DISCRETE FOURIER TRANSFORM 

In view of the interest of the bases {\ja; 0a) : a = 0, 1, . . . , 2j} for quantum information 
and quantum computation, we shall continue with the case r = 0. From now on, we shall 
also assume that a — 0, 1, . . . , 2j. Furthermore, by making the following change of notation 



n = j + m, \n) = |j, m), d = 2j + 1, 



(12) 



Eq. ([9]) gives 



\ja; 0a) = ^J2 q< d - n)a/2+na \n) . 
yd n =0 



Alternatively, the change of notation 



(13) 



k = j-m, \k) = \j,m), d = 2j + l 



(14) 



leads to 



d-l 



yd k=0 



Equations ( 1131) and (1151) were used in 



23 



(15) 



241 ] . They are equivalent as far as 



221 ] and 

quadratic discrete Fourier transforms and mutually unbiased bases (MUBs) are concerned. 
Both Eqs. f|T3|) and f|T5|) correspond to quantum quadratic discrete Fourier transforms which 
can be inverted to give 



1 



d-i 



\n) = —q-< d - n ) a l 2 ]T q- an \ja; 0a) 
yd a=0 



and 



d-i 



\k) = J- q -(k+D(d-k-i)a/2 £ q ^+i)\ ja .Q a y 
yd a=0 



(16) 



(17) 



6 



Note that the word quantum in quantum quadratic discrete Fourier transform refers to the 
fact that the vectors \n) or \k) (corresponding to \j,m)) are used in the quantum theory of 
generalized angular momentum. 

In the following we shall adopt the change of notation ( Fl2|) and shall re-define \ja; Oa) as 
\aa). In other words 

i d-l 

\aa) = -= £ q n ^ a / 2+na \n) (18) 
yd n=o 

or, in an equivalent way, 

i d-l 

\n) = q -< d - n ) a / 2 J2 q~ an \aa). (19) 
yd a=0 

Note that the action of the operator v 0a on the vector \n) reads 

v 0a \n) = q- {n+1)a \n + 1} (20) 

modulo d. It is clear that the basis 

B 0a = {\aa):a = 0,l,...,d-1} (21) 

is an alternative to the basis = B 2 j +1 . There are d — 2j + 1 bases of this type for a in 
the ring Z/dZ. Each basis B^ a spans the regular representation of the cyclic group Cd- 

IV. QUADRATIC DISCRETE FOURIER TRANSFORM 

The expression 

( tt \ _ 1 „n(d-n)a/2+na fooN 
y H 0a) na - —f=Q , [22) 

occurring in ffTK]) . defines the na th matrix element of a quadratic discrete Fourier transform. 
To be more precise, for fixed d and a, let us consider the transformation 

x = {x(n) eC:n = 0,l,...,d-l}<H-2/ = {y(a) e C : a = 0, 1, . . . , d - 1} (23) 

defined by 

d-l d-l 



y(a) = J2 (H 0a ) na x(n) & x(n) = ^ (H 0a ) na y(a). (24) 



n=0 a=0 
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The particular case a = corresponds to the ordinary discrete Fourier transform which 
satisfies 



(H 00 ) A = I d , 



(25) 



where Id is the identity d x d matrix. For a ^ 0, the bijective transformation x -H- y can be 
thought of as a quadratic discrete Fourier transform. The analog of the Parseval-Plancherel 
theorem for the usual Fourier transform can be expressed in the following way. 

Result 2. The quadratic discrete Fourier transforms x -B- y and x' <->■ y' associated with 
the same matrix H 0a , a G Z/dZ, satisfy the conservation rule 



d-i 



d-i 



V( a )y'( a ) = S x{n)x\n), 

a=0 n=0 



(26) 



where the common value is independent of a. 



It is to be observed that the matrix H 0a is a generalized Hadamard matrix in the sense 
that the modulus of each of its matrix element is equal to 1/ \fd. Such a matrix reduces the 
endomorphism associated with the operator t> 0a - As a matter of fact, we have 

(q 1 ... 0\ 
q 2 ... 



(H J V 0a H 0a = q^ 1 ^ 2 



V o o 







q d J 



(27) 



where the matrix 



Vr 



Oa 



{0 q a 
q 2a 





\ 


Jd-l)o 



(28) 



\1 ... / 

represents the linear operator v 0a on the basis 



B d = {\n) :n = d-l,d-2,...,0}, 



(29) 



known as the computational basis in quantum information and quantum computation. 
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The Hadamard matrices H 0a and (a,b G Z/dZ) are connected to the inner product 
(aa\bf3). In fact, we have 

-I d-1 



((H 0a ^H ob ) a 



(aa\b(3) = - g^-«)( 6 - a )/ 2 +"^- Q ). 



(30) 



n=0 



Thus, each matrix element of (Ho n )^ H nh can be written in the form of a generalized quadratic 
Gauss sum S(u, v, w) defined by |25j ] 



S(u,v,w) = e 



iir(un +vn)/w 



(31) 



n=0 



where u, v, and w are integers such that u and w are mutually prime, uw ^ 0, and uw + v 
is even. In detail, we obtain 

1 



(aa\b/3) = ((H 0a ^H ob ) = ^-S(u,v,w) 



with the parameters 



u = a — b, v = —(a — b)d — 2(« — f3), w = d, 



which ensure that uw + v is even. 

The matrix V 0a can be decomposed as 



Voa — X Z a , 



(32) 



(33) 



(34) 



where 



and 



/0 1 

1 



Vi o o 



(1 

q 

q 2 

\0 





1 

0/ 

o \ 






q d-lJ 



(35) 



(36) 
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The unitary matrices X and Z g-commute in the sense that 

X Z - qZX = 0. (37) 

In addition, they satisfy 

(X ) d = Z d = I d . (38) 

Equations (1371) and (I3"8"j) show that X (to be noted as X in what follows in order to conform 
to the notations used for Pauli matrices) and Z constitute a Weyl pair. Weyl pairs were 



introduced at the beginning of quantum mechanrcs 26] and nsed for building operator 
unitary bases [12]. The Weyl pair (X, Z) turns out to be an integrity basis for generating a 
set {X a Z : a,b G Z/dZ} of d 2 generalized Pauli matrices in d dimensions (see for instance 



23, 



24j and references therein). In this respect, note that for d = 2 we have 

X = a x , Z = a z , XZ = -ia y , X°Z° = a , (39) 

in terms of the ordinary Pauli matrices a = J 2 , a x , <r y , and a z . Equations (!37|) and ( 138|) 
can be generalized through 

V 0a Z - qZVoa = 0, ( V 0a ) d = e^ d - 1)a I d , Z d = I d , (40) 

so that other pairs of Weyl can be obtained from Vo a and Z. 

V. MUTUALLY UNBIASED BASES 

From a very general point of view, let us recall that two orthonormal bases B a = {\aa) : 
a — 0, 1, . . . , d— 1} and B\, = {\b(3) : /3 — 0, 1, . . . , d— 1} of the Hilbert space C d are said to 
be mutually unbiased if and only if the inner product (aa\b(3) has a modulus independent 
of a and (3. In other words 

Va e Z/dZ, V/3 G Z/ciZ : |(aa|6/3)| = r5 a , 6 r5 Qi/3 + (1 - 5 a , 6 )-^. (41) 

From Eq. fj4T|) . note that if two MUBs undergo the same unitary or antiunitary transfor- 
mation, they remain mutually unbiased. It is well-known that the maximum number M of 
MUBs in C d is M = 1 + d and that this number is attained when d is a prime number p 
or a power p e (e > 2) of a prime number p \\^-[la^. In the other cases (d ^ p e , p prime 
and e integer with e > 1), the number M is not known although it can be shown that 
3 < N < 1 + d (see for example [3]). In the general composite case d = riiPfS we know 



that 1 + min(p^) < Af < 1 + d (see for example [17|) 
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A. Mutually unbiased bases for d prime 

For d — 2, it can be checked that the bases B 00 , B i, and B 2 are 1 + d = 3 MUBs. A 
similar result follows for d = 3: the bases Bqq, Bqi, Bq2, and B% are 1 + d — 4 MUBs. This 
can be generalized by the following result. 

Result 3. For d = p, with p a prime number, the bases B 00 , B mi . . . , -B 0p _i, B p form a 
complete set of 1 +p MUBs. The p 2 vectors \aa) , with a, a — 0,1, ... ,p — 1, of the bases 
B 00 , B i, . . . , £>o P -i are given by a single formula (namely Eq. (TT^j j. 

The proof of Result 3 is as follows. First, Eq. ( {TBI) yields 

\(k\aa)\ = ^-, (42) 

a relation that holds for all k, a, and a in Z/pZ so that each basis is unbiased with B p . 
Second, the generalized quadratic Gauss sum S(u,v,w) in ( |32l) . with d = p prime, can be 
calculated to give 

|( aa | &/3 )| = _L (43) 

for all a, 6, a, and /3 in Z/pZ. This completes the proof. 

We note in passing that, in the case where d = p is a prime integer, the product (i^oa)^ -f^oft 
is another generalized Hadamard matrix. 

To close this subsection, we may ask what becomes Result 3 when the prime integer p 
is replaced by an arbitrary (not prime) integer d. In this case, the formula (|T8l) does not 



provide a complete set of 1 + d MUBs. However, it is possible to show [23|, |24j that the bases 
Boa, B 0a& i, and B d are three MUBs in C d (the addition © is understood modulo d). This 
result is in agreement with the well-known result accordin g to which the maximum number 
of MUBs in C d , with d arbitrary, is greater or equal to 3 ( la]). Moreover, it can be proved 



241 ] that the bases B 0a and B 0a92 are unbiased for d odd with d > 3 (d prime or nor 



prime) . 



B. Mutually unbiased bases for d power of a prime 



Equation (|T8|) can be used for deriving a complete set of 1 + p e MUBs in the case where 
d = p e is a power (e > 2) of a prime integer p. The general case is very much involved. 
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Hence, we shall proceed with the example p = e = 2 corresponding to two qubits. 

For d = 2 2 = 4, the application of ([TBI and f[2TT) yields four bases B 0a (a = 0, 1, 2, 3). As 
a point of fact, the bases -Both Bqi, Bq2, B03, and -B4 do not form a complete set of 1 + d — 5 
MUBs. However, it is possible to construct a set of five MUBs from repeated application of 

dUD. 

Four of the five MUBs for d = 4 can be constructed from the direct products \aa) © \b/3) 
which are eigenvectors of the operators v 0a <S> fob- Obviously, the set 

B 0a0b = {\aa) ® \b(3) : a, = 0,1} (44) 

is an orthonormal basis in C 4 . It is evident that -Boooo and -B0101 are two unbiased bases 
since the modulus of the inner product of \0a) © 1 0/3) by |la') <g> |l/3') is 

|(Oa|k*')(0/3|l/3')| = ^. (45) 

A similar result holds for the two bases -B0001 an d -B ioo- However, the four bases -Booooj 
-B0101, -B0001, and -B0100 are not mutually unbiased. A possible way to overcome this no-go 
result is to keep the bases -Boooo and -B0101 intact and to re-organize the vectors inside the 
bases -B0001 and -B0100 in order to obtain four MUBs. We are thus left with four bases 

Woo = Boooo, W u = B 0W i, W Q1 , W 10 , (46) 
which together with the computational basis B4 give five MUBs. In detail, we have 

Woo = {\0a)®\0P):a,P = 0,l}, (47) 

W n = {|la>® 1 1/3) :a,/3 = 0,l}, (48) 

W01 = {\\0a}® \ip)+n\0a®l)® |1^©1) :a,/3 = 0,l}, (49) 

W10 = {A|la)® \0P)+n\la®l)® |0/3©1) :a,/3 = 0,l}, (50) 

where 

1 _ i 1 + i /vis 

A = — , » = — ( 51 ) 

and the vectors of type \aa) are given by the master formula (TlBl . As a resume, only two 
formulas are necessary for obtaining the d 2 = 16 vectors \ab; af3) for the bases W a b, namely 

M^ooj^h : |aa;a/3) = \aa) <S> |a/5) , (52) 
W 01 , W w : |aa © 1; a/3) = \\aa) ® \a © 1/3) + /x|aa © 1) © |a © 1/3 © 1), (53) 
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for all a, a, and in Z/2Z. 

It is to be noted that the vectors of the Wqq and Wn bases are not intricated (i.e., each 
vector is the direct product of two vectors) while the vectors of the Wqi and W\o bases are 
intricated (i.e., each vector is not the direct product of two vectors). 

Generalization of (1521) and (1531) can be obtained in more complicated situations (two 
qupits, three qubits, . . . ). The generalization of (|52|) is immediate. The generalization of ( |53|) 
can be achieved by taking linear combinations of vectors such that each linear combination 
is made of vectors corresponding to the same eigenvalue of the relevant tensor product of 
operators of type vo a . 



VI. MUTUALLY UNBIASED BASES AND UNITARY GROUPS 



In the case where d is a prime integer or a power of a prime integer, it is known that the 
set {X a Z b : a,b = 0, 1, . . . , d — 1} of cardinality^ 2 can be partitioned into 1 + d subsets 



containing each d — 1 commuting matrices (cf. 15[). By way of illustration, for d = p with 



p prime, the 1 + p sets of p — 1 commuting matrices are easily seen to be 

V = {X°Z a :a = l,2,...,p-l}, (54) 

Vi = {X a Z°:a = l,2,...,p-l}, (55) 

V 2 = {X a Z a :a = l,2,...,p-l}, (56) 

V 3 = {X a Z 2a :a = l,2,...,p-l}, (57) 

V p -i = {X a Z^ a :a = l,2,...,p-l}, (58) 
V p = {X a Z ip - 1)a : a = 1,2,..., p- 1}. (59) 

Each of the 1 + p sets Vo, V±, . . . , V p can be put in a one-to-one correspondance with one 
basis of the complete set of 1 + p MUBs. In fact, Vo is associated with the computational 
basis while Vi, V2, . . . ,V P are associated with the p remaining MUBs in view of 

V 0a eV am , a = 0,l,...,p-l. (60) 

Keeping into account the fact that the set {X a Z b : a, b = 0, 1, . . . ,p — 1} \ {X°Z } spans 
the Lie algebra of SU(p), we get the following result. 
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Result 4. For d = p, with p a prime integer, the Lie algebra su(p) of the group SU(p) 
can be decomposed into a direct sum of 1 + p abelian subalgebras each of dimension p — 1, 
i.e. 

su(p) ~ v t±J vx W ■ ■ ■ tfcl v p (61) 

where the 1+p subalgebras Vo, v x, ■ ■ ■ , v p are Cartan subalgebras generated respectively by the 
sets Vo, Vx, ■ ■ ■ ,V P containing each p — 1 commuting matrices. 

Result 4 can be extended when d = p e with p a prime integer and e an integer (e > 2): 
;here exists a decomposition of su(p e ) into 1 + p e abelian subalgebras of dimension p e — 1 



27(-|33| (see also [23j,|24j) 



jmposi 
23|, I24 



VII. CLOSING REMARKS 

MUBs prove to be useful in classical information theory (network communication proto- 
cols), in quantum information theory (quantum state tomography and quantum cryptogra- 
phy) , and in the theory of quantum mechanics as for the solution of the Mean King problem 



241 ] for an extensive 



and the understanding of the Feynman path integral formalism (see [23j, 
list of references). 

There exist numerous ways of constructing sets of MUBs. Most of them are based on 
discrete Fourier transform over Galois fields and Galois rings, discrete Wigner distribution, 
generalized Pauli operators, mutually orthogonal Latin squares, discrete geometry methods, 
angular momentum theory and Lie-like approaches. In many of the papers dealing with 
the construction of MUBs for d a prime integer or a power of a prime integer, the explicit 
derivation of the bases requires the diagonalisation of a set of matrices. 

In the present paper, the generic formula (fT8l) arises from the diagonalisation of a single 
matrix (the matrix Vo a ), for the M = 1 + p MUBs corresponding to d = p with p a prime 
integer. Repeated application (e times) of this formula can be used in the case where d = p e 
is the power of a prime integer. Results 1 and 3 of this paper concern the closed form formula 
( |T8l) . Its derivation is based on the master matrix Vo a . From Vo a we can deduce the Weyl 
pair (X, Z) through 

X = V o, Z=(V m ) ] V m . (62) 
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The operators X and Z are known as the flip or shift and clock operators, respectively. For 
d arbitrary, they are at the root of the Pauli group, a finite subgroup of order d 3 of the group 
SU ( d) , of considerable importance in quantum information and quantum computation (e.g., 
see 34|- 36| and references therein for recent geometrical approaches to the Pauli group). 
The Pauli group is relevant for describing quantum errors and quantum fault tolerance in 
quantum computation. 
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